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Abstract: 

The ground and first excited state of two superconducting layers in interac- 
tion are studied considering two different coupling terms, one represented by 
the standard Josephson interaction, and one new, which is a superexchange 
pairing force between bilayer pairs. It is shown that a moderate-to strong 
Josephson interaction produces a low-lying collective state, pictured as an 
out-of-phase oscillation of the BCS gauge angles of the two layers. This anti- 
symmetric angular oscillation might explain the 41 meV resonance observed 
in the neutron scattering experiments. The bilayer pairs are formed by elec- 
trons from different layers with an antiparallel orientation of the spins, being 
related to the antiferromagnetic arrangement. The pair operators within 
the layers together with the bilayer pairs generate by commutation an so(5) 
algebra. It is shown that the transition between the superconducting and 
antiferromagnetic phases can be explained assuming the dependence on con- 
centration of the bilayer pairing strength, with maximum at half-filling. 



PACS numbers: 71.10.Li, 74.50.+r, 74.72.-h, 03.65.Fd 
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I. Introduction 



Since the discovery of the high-T c superconductivity ffl the remarkable 
properties of the cuprate compounds have been a subject of thorough experi- 
mental and theoretical investigation. Although a generally accepted theory is 
not yet established, the extensive knowledges accumulated during the years 
indicate that there are few key elements closely related to the mechanism of 
the high-T c superconductivity 0. 

A common element in all cuprate compounds is the layered structure, 
consisting of CuO planes separated by block layers ||. There is also a single 
parameter, the occupation fraction, whose variation may switch the state 
between superconducting (SC) and antiferromagnetic (AF). For example, 
when the planar Cu atoms of La2Cu04 or YBa2Cu30g are divalent each site 
contains one particle (half-filling), and the compounds are AF insulators. 
However, if an increasing number of Cu atoms become trivalent, (hole dop- 
ing), the AF order is destroyed, and after a transitory spin-glass state the 
compounds becomes metallic, or below T c , superconducting ||. 

Although the pair binding energy is supposed to come mainly from the 
interactions within the planes, the interlayer coupling is particularly impor- 
tant. Thus, the properties of the high-T c compounds become more isotropic 
in the SC state pf, and the measurements indicate that the current normal 
to the planes is of Josephson type ||. In the insulating phase, the AF order- 
ing is not restricted to the nearest neighbors within the planes. It appears 
also between the adjacent layers || 0, and is due to a weak AF interlayer 
exchange interaction [H, [J. 

An interlayer coupling of Josephson type can amplify the pairing corre- 
lation within the CuO planes, explaining the high T c -values and the con- 
densation energy |T0|, [11]. With this coupling the layered cuprates resemble 
a stack of Josephson junctions p2[, each of them being represented by two 
adjacent layers (bilayer) coupledoy the Josephson interaction. 

The low- lying spectrum of the bilayer system is of particular interest, and 
when the layers become superfluid one can expect collective excitations re- 
lated to the phase of the BCS (Bardeen-Cooper-Schreiffer) order parameter 



13|, pj| . As it was shown for the first time in JL4|, a time- dependent mean- 
field treatment (TDMF) of the Josephson coupling in a two-component SC 
system indicates the occurrence of a low-lying resonance. This is antisym- 
metric between the two components, and can be pictured as an oscillation 
of the relative phase angle between the two BCS condensates. However, 
the TDMF approximation is quasi-classical, and therefore in this approach 
purely quantum effects as tunneling and quantum coherence oscillations are 
usually suppressed |15| . 

The relevance of the TDMF phase angle resonance for the low-lying spec- 
trum of a SC bilayer system will be studied in Sect. I. Each component of the 
bilayer will be described by a reduced BCS interaction term, and spectrum is 
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calculated numerically as a function of the strength of the Josephson coupling 
and the occupation fraction. It will be shown that for moderate-to-strong 
coupling the phase angle resonance gives a very good description of the first 
excited state. This resonance is antisymmetric between the planes, appears 
only in the superconducting state, is coherent between the CuO planes and 
of nonmagnetic origin. By these properties it can provide an explanation 
for the 41 meV resonance observed in the neutron scattering experiments on 
YBa 2 Cu 3 7 MM- 

If the two layers can exchange particles, then a pair of electrons localized 
in one plane may be changed into a pair of electrons from different planes. 
Such "bilayer pairs" appear also when the Josephson coupling term between 
the two layers is derived from a general two-body pairing interaction. Be- 
cause these pairs contain electrons from different planes with an antiparallel 
arrangement of the spins, they can be naturally related to the AF inter- 
layer ordering. Following the example presented in |14|, the extension of the 
su(2)©su(2) algebra of the BCS pairing operators within the planes by the 
bilayer pair operators leads to the so(5) algebra, and it will be presented in 
Sect. III. 

A model of the high-T c superconductors where so (5) is generated by the 
three operators of the total spin, the particle number operator, and six pair 
operators which generate "rotations" between the total SC order parameter 
and the Neel vector was formulated in |18], |I9|| . The AF state at half-filling 
appears as a solid formed by Cooper pairs, and AF-SC transition corre- 
sponds to the melting of this solid under the influence of a fictious "B-field", 
represented by the chemical potential However, in this approach the 

separation of the system in the superconducting components associated to 
the CuO planes is not explicit, and it is difficult to establish the connection 
between the model parameters and the strength of the interlayer Josephson 
coupling. Therefore, the 41 meV excitation observed in YBa2Cus07 appears 
as a Goldstone mode associated to the breaking of the SO (5) symmetry rather 
than to the "restoring force" determined by the Josephson interaction, and 
its generators have both SC and AF components [pOfl . 

The pairing interaction between the bilayer pairs considered here resem- 
ble the superexchange interlayer interaction 0, and is complementary to the 
Josephson coupling. The properties of the low-lying spectrum of the bilayer 
system will be studied as a function of the strength of this interaction and 
the occupation fraction in Sect. IV. It will be shown that by contrast to the 
Josephson coupling, the new interaction may suppress the superconductivity 
of the planes, leading to an AF arrangement. A summary of this study and 
the conclusions are presented in Sect. V. 
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II. Low-lying excitations induced by the Josephson coupling 

The model Hamiltonian considered here is supposed to describe the elec- 
trons of two close CuO layers (bilayer system), denoted A and B. Assuming 
that these layers are related by a translation along the vector Q, between 
their lattice sites there is a one-to-one correspondence expressed by b = a+Q. 
The electrons of each CuO plane are supposed to interact strongly by pairing 
forces, and within a restricted BCS approximation, the Hamiltonian has the 
form 

H 0L = eN L - G L PtP L , L = A, B. (1) 

Here 

N A = E( C 1t C «T + C li C ai) , N B = E( C 1+QT C «+QT + C i+Ql C a+Ql) , (2) 
a a 

are the particle number operators, 

p1 = E - «')i T i'i > p l = T,H^-^)ci +Q ^ a/+Ql , (3) 

a,a' a, a' 

create electron pairs with the wave function $(a — a') = $(a' — a), and c£ m 
(cam) are the creation (annihilation) operators for an electron at the site a 
(a + Q) of the lattice A (B) with the Z-component of the spin up (m =|) or 
down (m =|). The layers are coupled by a Josephson interaction term 

Hj = -Gj{P ] a P b + P ] b P a ) , (4) 

such that the total Hamiltonian has the form 

H = e A N A + e B N B - G A P\P A - G B P B P B - Gj(P\P B + P B P A ) . (5) 

The model parameters are the single-particle energies e A , e B , the strengths 
G A , G B , Gj of the pairing forces and the number of the lattice sites. 

The commutation relations between the pair creation operators and their 
Hermitian conjugates are determined by the orbital wave function $(a — a'). 
For a d x 2_ y 2 -pair field, this function can be expressed as |4j] 

$( a - a) = ~{5 a -a',x + ^a-a'-x ~ 8a-a',y - ^a-a'-y) (6) 

where x, y are vectors along the x and y axes of the plane lattice relating the 
nearest neighbors. The commutator [P A , P A ] contains the sum 

$(a - ai)$(aa - a') = 5 a ^ + -(5 a ^ a > i2x + § a -a'-2x + $a-a>,2 y (7) 

ai 4 
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~^~&a—a',—2y 25 a — a i ,x+y 25 a — a ' t — x —y 25 a — a '^ x —y 2i3a—a',—x+y) 

Here only the first term of the right-hand side, which is diagonal, relates 
nearest neighbors. Assuming that the contribution due to the other terms 
can be neglected, and that the number of lattice sites Q is the same for A 
and B, then [P[, P L \ = N L - SI, L = A,B. Defining P 0L = (N L - Vt)/2, 

then [Pol, P}] = Pli an d eacn se ^ °f three operators (P[, Pql, Pl), L = A,B, 
generates an su(2) algebra. 

If Gj = the ground state \g(N p ) > of H is non-degenerate when the 
total number of pairs N p is even, and double degenerate when it is odd. 
Thus, \g(2k) > can be factorized in a product of two A;-pair states, one for 
each layer, as \g(2k) >= \A(k) > \B(k) >, where 

\L(k) >= )j J^^(Pt) k \° > ' L = A,B (8) 

and |0 > is the particle vacuum. When N p = 2k + 1 the ground-state 
doublet states are \g a (2k + 1) >= \A(k + l)B(k) > and \g b (2k + 1) >= 
\A(k)B(k + 1) >. The term Hj rise the degeneracy between these states, 
leading to the entanglement of the two planes. When Gj is small, the ground 
and first excited state of the coupled system are 

\g w (2k + l) >= -±={\A(k + l)B(k) > +\A(k)B(k + l) >) , (9) 
v2 

and 

\l w >=^=(\A(k + l)B(k)>-\A(k)B(k + l)>) . (10) 
v2 

The difference between their energies E x = E 1 — E is given in this case 
by the tunnel splitting E t = 2\ < A(k + l)B(k)\Hj\A{k)B(k + 1) > | = 
2GjVl 2 (n + Q,/2)(l — n + Sl/2) , where n = N p /2Vt denotes the occupation 
fraction. 

This result was compared with the excitation energy of the first level 
calculated numerically for a bilayer system with Q = 12, €a = £b = Q and 
Ga — Gb — 1, as a function of n and Gj. The first excited state is always 
antisymmetric between the planes, and when Gj is smaller than ~ 0.05, E x 
shows a strong odd-even effect at the variation of N p . This behavior is il- 
lustrated for Gj = 0.005 in Fig. 1 (a) (dots joined by dashed lines). When 
Gj > 0.1, this strong odd-even effect disappears, as it is shown in Fig. 1 (c) 
for Gj = 0.2. The comparison with E t (Fig. 1 (a,c), solid line) shows a good 
agreement when N p is odd and Gj is small, but for N p even or when Gj 
increase the differences become large, especially at half-filling. This result 
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indicates that at moderate-to-strong coupling (Gj > 0.1) the tunneling ap- 
proach to the Josephson interaction is not realistic. It is interesting to note 
that the range of the Gj values when the coupling is "weak" or "strong" 
depends of Q. For instance, decreasing Q from 12 to 6, the interval sepa- 
rating the weak and strong coupling regimes changes from ~ (0.01,0.1) to 
-(0.05,0.2). 

The pair transfer matrix element between the ground states, A = G < 

g(Np)\Pl\g(N p — 1) > is a measure of the "softness" of the system with re- 
spect to the variation of the number of pairs. A high value of A indicates the 
increase of the two-particle correlations in the ground state, and may signal 
the transition to a superconducting state. This parameter was calculated for 
Gj = 0.005 and 0.2, and the results are represented in Fig. 1 (b) and (d), 
respectively. Their comparison confirms the increase of the pairing correla- 
tions by the interlayer Josephson interaction, and shows that this effect is 
stronger at half-filling. 

When a layer becomes superconducting, its ground state breaks the U(l) 
symmetry of the Hamiltonian generated by the particle-number (charge) op- 
erator, and is well approximated by the BCS function 

\BCS(<f>,\) >=e (2 ^- 2z * Pi) |0> ,L = A,B . (11) 



Here A is the Fermi energy, <ft is the BCS "gauge angle", z = pe and p 

is a variable related to the BCS gap A = G < BCS(0, X)\Pl\BCS(0, A) > 
by A = (GQ/2) sin4p, and to the occupation fraction n =< BCS(0, A)| 

N L \BCS(0, A) > /2Q by n = sin 2 2p. The spontaneous breaking of the 
U(l) symmetry generates a specific collective phase-space |21| parameterized 



by the BCS "gauge" angle cf> and the canonically conjugate momentum p 
Qn = Q/2 — (e — X)/G ||14j| . In the absence of the interlayer coupling, <p has a 
linear dependence on time, 4>(t) = 2Xt, such that the system performs a free 
rotation in the BCS gauge-space, with the angular velocity 2A. This rotation 
can be interpreted as the Goldstone mode related to the U(l) symmetry 
breaking, and the classical Hamilton function associated to this evolution is 



defined by the expectation value of H l, HI 



H ((P,p)=<BCS\H 0L \BCS>=2(e-G^)p + Gp 2 . (12) 



This expression indicates that the inertial parameter corresponding to the 
free gauge rotation is Iq = 2/G. 

When the two superfluid layers are coupled by the Josephson interaction 
term Hj, the Hamilton function becomes 

H{<f>A,PA,<f>B,PB) =< BCS A BCS B \H 0A + H 0B + Hj\BCS A BCS B >= (13) 
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2(e-G^)(pA+PB)+G(p 2 A +p 2 B )-2Gjy/pAPB(£l - PaHV - Pb) cos(0 a -0 b ) . 

In this case, the dynamics of the system consists of a superposition between 
a uniform rotation of the "center of mass" angle {<pA + 4>b)/2, symmetric 
between the planes, and an oscillation in the relative angle (4>a — 4>b)/2, 
which is antisymmetric. Assuming the same occupation fraction n for both 
planes, in the harmonic approximation the excitation energy of this mode is 



E r = V2kC . (14) 

Here k = 2{G — Gj) + Gj/n(l — n) is the related inertial parameter, and 
C = 2Gj(A /G) 2 = 2GjVL 2 n{\ - n) « E t , is the strength of the restoring 
force determined by the Josephson coupling. The same excitation energy is 
obtained when Hj is treated within the random-phase approximation (RPA) 
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.He energy E r of the antisymmetric resonance oscillation and the BCS 
gap Ao obtained with the same parameters as used in the numerical calcu- 
lations are represented by solid lines in Fig. 1 (a,c) and (b,d), respectively. 
The agreement with the numerical values is very good at moderate to-strong 
coupling, when E x Pd E r and AkA . At weak coupling E r still reproduces 
the average trend of E x , but A is systematically smaller than A by a fac- 
tor ~ 0.83. Both A and A measure the pair correlation functions within 
a single plane. If each plane develops its superconductivity independently, 
then in the limit Gj — > one would expect that A — > Ao- However, this 
is not observed because when G j is very small, at N p odd the coupling be- 
tween the two planes leads invariably to an entangled ground state of the 
form \g w >. Therefore, the one-pair transfer matrix element A is reduced 

to l/y/2 of the value corresponding to an isolated plane. By contrast, at 
moderate-to-strong coupling the ground state is spread over many config- 
urations, and BCS approximation becomes realistic. This shows that the 
interlayer Josephson coupling has an important contribution to the occur- 
rence of the superconductivity within the planes. 

Beside the "kinetic" term G{p 2 A +p 2 B ) determined by the pairing forces, an 
additional contribution to 7i appears from the electrostatic energy between 
the planes, E es = e 2 (pA—PB) 2 d/{2e S) ( 0, p. 124), where d is the interlayer 

spacing and S the layer surface. With this term E r becomes \j2(k + k e )C, 

where k e = 2e 2 d/ (eoS). However, by contrast to k, which is size-independent, 
k e vanishes at the macroscopic limit S — > oo, and it will be neglected. 

The condensation energy and the superconducting energy gap are known 



from the measurements of the electronic specific heat [25J . For YBa2CusO 
the zero temperature superconducting gap has the maximum A = 20 meV, 
and the condensation energy is U c = 0.5 meV/f.u. cell. Assuming that 



A = GQ,Jn{l — n), with n ~ 1/2, and Gj/G w 0.5, one obtains E r 50 
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meV. This low value suggests that the antisymmetric angular oscillations in 
the BCS gauge space presented here may provide an explanation for the 41 
meV resonance observed in the neutron scattering experiments TO , |T7j| . 

The general form of the operator creating a pair of electrons with a three- 
dimensional orbital wave function ^(7 — I') in the bilayer system is 

p^ = Y,m-i')4Ai , (15) 

i,v 

where / runs over the lattice sites of both planes, A and B. Assuming that 
the in-plane section \l/(a — a') is the same as $(a — a'), and ^(7) is periodic 
in the Q-direction, \I/(Z) = ^{l + Q)^ this general pair creation operator can 
be decomposed as 

pt = pt + pt + pt ( 16 ) 

where P| and Pg are the pair creation operators within the planes of Eq. 
(3), and 

fit = £ *(a' - a)(ci +0T 4, i + ci +Qi cl, T ) (17) 

a, a' 

creates a pair of electrons localized in different planes, with an antiparallel 
spin orientation. 

The antiparallel orientation of the nearest spins from different planes is 
a characteristic of the inter layer antiferromagnetic ordering. In fact, the AF 
ordering of the high-T c materials was described within the resonance valence 



bond (RVB) theory as a condensate of nearest neighbors singlet pairs [p6 
A general pairing interaction of the form —GP^P leads naturally the terms 
of the Hamiltonian H, but contains also other terms dues to the "bilayer 
pairing" operators R. In particular, a pairing interaction —R^R between the 
bilayer pairs can represent a convenient approximation for the superexchange 
interlayer coupling presented in [2[, p. 395. This term can produce the AF 
condensation, and it will be studied by extending to so(5) the su(2)©su(2) 
algebra used to describe the effect of the interlayer Josephson coupling. 

III. The so(5) algebra of the pair operators in the bilayer system 

The shift along the chain 

pt _» pt _ pt (18) 

of the three pair creation operators defined above is generated by commuta- 
tion with the interlayer "hopping" operator 

T + = E( c i+Ql c <4 - C 1+QT C <4) > ( 19 ) 
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such that 

[T + ,P\] = -Rj , [T +) tf]=2P B , [T + ,Pl] = . (20) 
Similarly, the down shift 

P t _ jet _> pt (21) 

is generated by commutation with T_ = Each set of three operators, 

(R\Ro,R) and (T+,T ,T_), with R Q = P 0A + P 0B and T = (iV B - N A )/2 
generate a new su(2) algebra, while the set of the 10 operators T + , T_, P4, P4, 
PJ,R, P b ,Pb, Poa,Pob generate the so(5) algebra. 

It is important to note that by contrast to P A and P B which create spin 
singlet pairs, the pair created by R^ is the m = component of a spin triplet. 
A spin singlet pair of electrons from different planes is generated by 

St = £ $( a - a')(ci +QT cl, i - ci +Qi cl, T ) . (22) 



Pairs of this type have been used to construct the RVB vacuum in [ 27] . The 
shift along the chain 

P\ - St - P B (23) 
is generated by the commutation with the hopping operator 

r + = E( c 1+qi c <4 + C l+QT C at) > ( 24 ) 

a 

and r + , r_, P4, P A , ^t, 5, Pg, Pb, Pqa, Pqb generate also an so(5) algebra. 

In the defining representation, the so(5) generators are 5x5 matrices 
fpq = £pq — £-q-p, where the indices p and q take the values -2,-1,0,1,2, and 
e pq is the 5x5 matrix with the element on the row p and the column q equal 
to 1, and all the rest 0. Thus [e pq , e^] = 8 q ke p i — Si p ek q , and the commutation 
relations of the so (5) generators are 

[fpq, fkl] — fiqkfpl — Splfkq + ^p-kf-lq ~ &-qlfp-k ■ (25) 

The correspondence between the generators with negative roots of the defin- 
ing representation of so(5) and the two different realizations provided by the 
pairing operators, denoted so(5) r and so(5) s , is summarized in the Table 1. 
One should note that these two realizations are related simply by a change 
of sign for all the single-particle states with spin up in the layer A. This 
transformation can be written as cL — > — cL, and as long as it leaves the 
Hamiltonian invariant, the results are independent on the choice of so(5) r or 
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so(5) s for the explicit realization. 

The Cartan subalgebra of so(5) is generated by two elements, /22,/n, 
and in both realizations 

f22 = l(P0B + P A) , fll = \(P B-P0A) ■ (26) 

The algebra so (5) is semisimple of rank 2, and its irreducible representations 
can be labeled by the numbers [Aq, N 2 ], N\ > N 2 > of the associated Young 
diagram, or by the eigenvalues of the generators of the Cartan subalgebra 
for the highest weight state PBf . If this state is chosen to be the particle 



vacuum, |0 >, then similarly to the case of su(3) J29|, the states \klm > 

\klm>=(f 2 „ 1 ) k (f 12 ) l (f 02 ) m \0> , (27) 

with k + I + m < 2Q span the Hilbert space 7Yq, carrier of the symmetric 
irreducible representation [Q, 0] of so(5). 

The labeling of the states of this irreducible representation requires three 
indices. One of them is provided by the total number of particles 2(k + l+m), 
which takes values between and 4f2, and another by the eigenvalue k — loi 
T Q (= r ). The choice of the third index is a difficult long-standing problem, 
because the set of all possible states \klm > with a fixed N and k — I is 
non-orthogonal and overcomplete [BOl. For instance, there is a non- vanishing 
amplitude < kk0\002k >= {2k)\k\\Wj [Vt — k)\ for recoupling 2k bilayer pairs 
in two sets of k one-layer pairs. 

The basis of the [Cl, 0] irrep spaces can be obtained by using vector co- 
herent states methods, by projection from a highest weight state |HJ. Here 



the basis will be constructed numerically, restricting the overcomplete set of 
Eq. (27) by the Gramm-Schmidt procedure. As a test, the dimension of the 
spaces l~Cfi generated for 2 < Q < 12 was found to be the same as given by 
the analytical formula 

d(n,0) = i(2ft + 3)(J] + 2)(ft + l) . (28) 
6 

The action of the generators on the (non-orthonormal) states of Eq. (27) is 
expressed by ] 2 -\\klm >= \k + llm >, f 2 i\klm >= \kl + lm >, f 2 o\klm > = 
\klm + 1 >, and 

f-i 2 \klm >= -k(k + m-Q-l)\k-llm > +m(m- l)\kl + lm -2 > (29) 
fi 2 \klm >= -1(1 +m - Q - l)\kl - lm > +m(m - l)\k + llm - 2 > (30) 

Tfl — 1 

f 02 \klm >= kl\k-ll-lm + l > -2m(k + l-Q+ )\klm-l > . (31) 
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These formulas are useful to calculate the matrix elements of the general 
pairing Hamiltonian in the so(5) basis. 



IV. SC-AF transition induced by the bilayer pair interaction 

The extended Hamiltonian which includes the interaction between the 
bilayer pairs has the general form 

H g = -G(P A P A + P B P B ) - Gj(P\P B + P B P A ) - grfR . (32) 

This Hamiltonian is expressed only in terms of the so(5) r generators, and 
therefore it allows to restrict the Hilbert space to a space Ti^ of irreducible 
representation. 

The low-lying spectrum of H g has been calculated for G = 1, Gj = 0.5 
and Q = 6, as a function of the strength g of the bilayer pairing interaction 
and the occupation fraction n. The ground state energy E and the exci- 
tation energy of the first level E x = E\ — E are represented as a function 
of n and g in Fig. 2(a) and 2(b). The pair transfer matrix element A is 
represented in Fig. 3(a). 

For the study of the AF correlations each lattice L = A, B is separated 
in two sublattices, denoted by (L, 1) and (L, 2), such that the nearest neigh- 
bors of the spins from (L, 1) within the plane L are in (L, 2). If the nearest 
neighbors of (A, 1), (A, 2) outside the A-plane are in (B, 1), {B,2), then in 
a classical antiferromagnetic configuration the spins from (A, 1) and (B, 2) 
have all the same orientation, opposite to that of the spins from (A, 2) and 

In a quantum system the situation is more complex, because the ground 
state is a symmetric superposition of classical AF states with opposite ori- 
entations of the Neel vector | 3~2|| . Thus, as a measure of the AF correlations 
in the quantum bilayer system here will be considered the ground state ex- 
pectation value of the operator C = (5^(1) - Sf{2)){Sf (1) - Sf (2)), where 

S t( a ) = \ £ (4t c «T - c l| c a|) , (33) 
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ad{A,a) 



and 

S B (ry) = 

2 



Sf(a) = ~ Wa+Q] C a+Q^ ~ c\ +Qi C a+Q[ ) , (34) 

ae(B,a) 



for a = 1, 2 are the z (c-axis)-components of the total spin for each sublattice. 
If only the terms containing nearest neighbors are retained, then the ground- 
state expectation value < g\C\g > can be approximated by < g\Sf{l)Sf (1) + 
Sf{2)Sf{2)\g > and also by < g\S?S?\g >, where = 5 2 L (l) + ^(2). This 
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approximation is particularly convenient if the Hilbert space is restricted to 
Tin- In this case, the matrix elements of S^Sf can be calculated analytically, 
and have the expression 

< k'l'm'\SfSf\klm >= [{k + j)(l + j)-kl- m k + l 2 +m + (35) 

m(m-l) (^ 2 -l) m(m- 1) 

2 fi(2fi-l)-l^ + n(2f2-l)-l P ' 

where / =< k'l'm'\klm > and g =< k'l'm'\k + 1/ + lm — 2 >. 

The expectation value of this operator in the ground state of H g is repre- 
sented as a function of n and g in Fig. 3(b). The results indicate that when 
g is small A > 0, < g\S^Sf\g >= 0, and the main components of the ground 
state are the states \klm > with k + l > m. However, when g increases above 
~ 0.8 the structure of the ground state has a sharp transition, and the states 
\klm > with m > k + l become dominant. During this transition A vanishes 
(Fig. 3 (a)) and < g\S^Sf\g > becomes negative (Fig. 3 (b)), indicating 
the onset of the AF ground-state correlations. In the transition region the 
excitation energy E x has maxima when the number of pairs is even, and min- 
ima when it is odd (Fig. 2(b)). This behavior is correlated with the change 
in the structure of the first excited state, whose main components can be 
obtained by turning symmetrically one bilayer pair of the ground state into 
a pair within the layers. 



IV. Summary and conclusions 



The interlayer coupling in the cuprate compounds has a fundamental role 
in the physics of the high-T c superconductivity, and in the recent years an 
increasing experimental and theoretical effort was devoted for understanding 
its mechanism 0. In this work the low- lying properties of two supercon- 
ducting layers in interaction have been described considering two different 
coupling terms, one represented by the standard Josephson interaction, and 
one new, which is a " superexchange" pairing force between bilayer pairs. 

The Josephson coupling between two superconductors is usually presented 
as a perturbation allowing the flow of the Josephson tunneling current []33| 



Complementary to this result, a TDMF |14[ or RPA ||22j| calculation indicates 
that with this interaction the difference between the BUS phase angles of the 
two superconductors is not a constant, but oscillates in time, producing an 
excitation which could be observed as an antisymmetric low-lying resonance. 
These two pictures have been compared (Fig. 1) with the exact results ob- 
tained numerically for a restricted BCS Hamiltonian. 

It was shown that the tunneling approach describes very well the energy 
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of the first excited state when the number of pairs is odd and the strength 
of the Josephson coupling is small. In this case the bare ground state is 
double degenerate. The Josephson interaction rises this degeneracy, split- 
ting the doublet into a non-degenerate ground state, symmetric between the 
two components, and the first excited state, which is antisymmetric. The 
two components of the bare ground state become non-stationary states, in- 
terchanged periodically in time by quantum coherence oscillations. These 
oscillations can be pictured as a periodic exchange of a pair between the two 
layers. One should note, however, that at the relatively high temperature of 
the cuprate superconductors they might be completely suppressed by ther- 
mal decoherence |34[ . 

At moderate-to-strong coupling the excitation energy of the first level is 
very well described as a collective out-of-phase oscillation of the BCS gauge 
angles of the two layers. By contrast to the phase of a quantum wave func- 
tion, the BCS angle is a collective variable defined only in the superconduct- 



ing state, canonically conjugate to the particle number [14], g9|. Assuming a 
ratio Gj/G = 0.5 between the strength of the Josephson coupling Gj and the 
pairing force within the layers G, the excitation energy of this resonance in 
YBa2Cu307 is ~ 50 meV. This is a promising result, because indicates that 
by a more accurate estimate of Gj/G, the antisymmetric angular oscillation 
in the BCS gauge space might explain the 41 meV resonance observed in the 
neutron scattering experiments JR], 0| . 

The effect of the Josephson coupling on the two-particle correlations was 
studied by calculating the pair transfer matrix element A between ground 
states. An increase of A indicates that the structure of the ground state is 
less sensitive to the addition or removal of a pair, and becomes closer to the 
one of the BCS condensate. The results presented in Fig. 1 (b) and (d) show 
that the Josephson coupling can induce the configuration mixing required for 
the occurrence of the superconductivity within the layers. At moderate-to- 
strong coupling A is very close to the BCS gap (Fig. 1 (d)), indicating that 
the ground state of the whole system can be well approximated by a product 
of two independent BCS functions, one for each plane. 

The decomposition of the most general pairing interaction in a bilayer 
system leads naturally to the Josephson coupling term, but also to a pairing 
interaction between "bilayer pairs", similar to the superexchange force dis- 
cussed in ref. ||. These pairs are formed by electrons from different planes 
with an antiparallel orientation of the spins, and therefore are related to the 
AF arrangement. The three different pair creation operators, (two for pairs 
within the layers and one for the bilayer pairs), together with their Hermi- 
tian conjugates generate by commutation the so(5) algebra. This result is 
important, because it shows that for a wide class of pairing Hamiltonians the 
many-body Hilbert space can be split in a series of finite-dimensional spaces 
of so(5) irreducible representation. In particular, one can expect to find the 
ground state within a space Hq, carrier of a symmetric irrep of so(5). The 
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construction of these spaces is a non-trivial problem in the Lie algebra rep- 
resentation theory, and it was presented in Sect. III. 

The low-lying properties of two superconducting layers coupled by Joseph- 
son and bilayer pair interactions have been studied as a function of the oc- 
cupation fraction in Sect. IV. For the explicit realization was used a set 
of so(5) r operators, but this choice is irrelevant as long as the Hamiltonian 
contains only a single type of bilayer pairs (R or S), and is invariant to the 
so(5) r <->so(5) s transformation. The results indicate that when the strength 
of the bilayer pair interaction g is ~ 0.8G or greater, there is a strong odd- 
even effect in the excitation energy of the first state (Fig. 2(b)), and the 
superconductivity within the planes disappears (Fig. 3 (a)). This transition 
is accompanied also by the onset of an AF arrangement between the spins of 
the two layers (Fig. 3 (b)), and is practically independent of the occupation 
fraction. Therefore, within this model an AF-SC transition at the variation 
of n can be explained only by assuming that the ratio g/G depends on n and 
reaches the maximum at half-filling. The increase of the pairing interaction 
strengths within the planes at half-filling was noticed in ref. ||35|| . The study 
of the concentration dependence of the interlayer coupling strengths appears 
therefore SIS db subject worth of further investigation. 
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Table 1. The so(5) generators with negative roots 
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T + /V2 


so(5) s 
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T+/V2 



Figure Captions 

Fig. 1. The excitation energy E x = E\ — Eq (a,c) and the pair transfer 
matrix element A = G < g(N p )\P\\g(N p — 1) > (b,d) as a function of the 
occupation fraction n = N P /2Q when Q = 12, Ga = Gb — 1, g — (dots 
joined by dashed lines). For comparison, E r , E t and Ao are represented by 
solid lines. Two values of the Josephson coupling strength are considered, 
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Gj = 0.005 (a),(b) and Gj = 0.2 (c),(d). 



Fig. 2. The ground state energy Eq (a) and the excitation energy 
E x = Ei — Eq (b) as a function of the occupation fraction n and the strength 
of the bilayer pairs interaction g when Q = 6, Ga = Gb = 1 and Gj = 0.5. 

Fig. 3. The pair transfer matrix element A (a) and the correlation 
function < g(N p ) S^S^\g(N p ) > (b) as a function of the occupation fraction 
n and the strength of the bilayer pairs interaction g when Q = 6, Ga = Gb = 
1 and Gj = 0.5. 
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